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ABSTRACT 


A laboratory study has been made of the characteristics 
of an interferometric type optical fiber pressure gradient 
hydrophone. The optical system is) configured as an 
all-fiber Mach-Zehnder interferometer excited by a single 
mode gallium arsenide diode laser. A pair of identical 
fiber coils, one in each arm of the interferometer, forms 
the sensing portion of the gradient hydrophone. Each coil 
consists of 10 meters of polyethylene jacketed, single mode 
fiber, wound in a doughnut shaped element of mean diameter 4 
cm and thickness 3 mm. An acoustic calibrator, similar in 
design to the USRD type G40 shipboard calibrator used by the 
U.S. Navy, was constructed for use in obtaining single and 
dual (dipole) coil sensitivities. The acoustic sensitivity 
of each coil was determined separately and then the pair was 
aligned coaxially and separated by 10 cm to operate as a 
pressure gradient device. Details of the construction of 
the system, calibration procedure, lock-in demodulator, and 
hydrophone sensitivity data are presented. The system has 
demonstrated a 14 dB sensitivity advantage over the USRD 
type G63 standard (piezoelectric) pressure gradient 


hydrophone. 
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I. INTRODUCTION 


A. BACKGROUND 


The use of optical fibers for communication and sensing 
is a rapidly expanding area of recent technology. Optical 
fiber technology is a little over a decade old, starting 
with the first demonstrations of low-loss (less than 1000 
dB/km) optical fiber waveguides. Technological advances in 
material processing, fabrication of optical fibers, 
interconnection devices, cables, sources, and detectors now 
allow repeaterless light transmission over great distances 
(tens of kilometers) with an extremely large bandwidth (Ref. 
ae 

Capitalizing on the availability of optical fiber 
components, there has been significant interest and progress 
in the past half decade in exploitation of optical fibers as 
transduction elements in anew class of sensors. These 
sensors are capable of detecting acoustic fields, linear and 
rotational acceleration, temperature, strain, electric and 
magnetic fields, and other physical parameters. This 
versatility can also lead to difficulties and much effort 
has been expended to enhance sensitivity to the parameter of 
interest while desensitizing the fiber to all other 


fluctuations CRef. 2J. Fiber sensors offer a number of 


il 


potential advantages: increased sensitivity, geometric 
versatility as fiber sensors can be configured in arbitrary 
shapes, lightweight components, and low cost. Progress in 
demonstrating these advantages has been substantial in the 
past few years with over 60 different types of sensors using 
optical fiber as the transduction element being developed. 
These sensors range from simple on/off fluid level 
indicators, to the more sophisticated interferometric 
configurations. This large number of individual devices is 
usually categorized into amplitude or phase sensors. In 
effect, the sensor modulates some parameter of the lightwave 
in an optical fiber, such as the intensity or the phase. “In 
the former case the physical perturbation interacts with the 
fiber to directly modulate the intensity of the light in the 
fiber. In the latter case the physical perturbation 
modulates the optical phase of coherent light propagating 
through the fiber. Optical interferometry must be utilized 
to convert the optical phase modulation to optical intensity 
modulation in order to permit photodetection. This 
technique allows for extremely high sensitivity. Phase 
modulation, conversion of phase modulation to intensity by 
interferometry, and subsequent detection by photodiodes will 


be discussed in more detail in Chapter II. 
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B. PURPOSE OF THIS STUDY 


Fiber optic acoustic sensors were first demonstrated to 
be feasible in 1977 (Ref. 3 and 4]. Since then a significant 
degree of progress has been achieved, both in the 
understanding of the relevant acousto-optic transduction 
mechanisms as well as in component development and sensor 
packaging CRef. 5 and 6]. Most of the fiber optic hydrophone 
wor k per formed to date has been with monopol ar 
omnidirectional devices. Figure 1.1 illustrates the basic 
geometry of atypical fiber optic hydrophone which utilizes 
acoustically induced phase modulation for transduction. a 
laser beam is split into two arms of an interferometer where 
a sensing coil in one arm is exposed to the acoustic 
pressure field while a reference coil, forming the second 
arm, is isolated from the field. The two beams are 
recombined and allowed to interfere on the surface of a 
photodetector. 

Since their invention in the late 1970‘s two unique 
characteristics of interferometric type fiber optic 
hydrophones, i.@., their basic two-arm design and their 
proven high sensitivity, have suggested their use as 
gradient-type devices. A fiber optic gradient hydrophone 
has a sensing coil in each arm of the interferometer = and 
both are exposed to the pressure field. Figure 1.2 


illustrates the basic geometry of an acoustic fiber optic 
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gradient hydrophone. It was the aim of this project to 
construct and partially evaluate in a laboratory environment 


the characteristics of a fiber optic gradient hydrophone. 
C. FORMAT OF THE REPORT 


The theoretical bases of conventional gradient 
hydrophones, the calibration of gradient hydrophones, and 
the behavior of interferometric sensors used in preparation 
of this report is presented in Chapter II. Details of the 
construction of the optical fiber sensor are presented in 
Chapter III. Included in Chapter III is a description of the 
characteristics of the individual parts of the 
interferometer. The parts are discussed in the order that 
they were built into the system. Specific instrumentation 
used to obtain data is also discussed in Chapter III. The 
experimental procedures used to establish the 
characteristics of the system and conduct sensitivity 
measurements of the individual sensor coils and the gradient 
(dual coil) hydrophone are discussed in Chapter IV. Analysis 
of the data and evaluation of the results are also presented 
in Chapter IV. Chapter V contains concluding remarks = and 
recommendations for further work. Appendix A lists a BASIC 
computer nguace program used in data acquisition and 
Appendix B presents an example of the data obtained using 


this program. 
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II. THEORY 


A. CONVENTIONAL GRADIENT HYDROPHONE 


In many acoustic applications it is necessary to 
determine the direction of signal arrival as well as to 
detect acoustic pressure levels. Directional features are 
usually achieved with the use of multiple hydrophones 
spatially distributed in some specified orientation, e.g., a 
line array. The simplest of these designs is the dipole 
pressure gradient hydrophone. 

Some pressure gradient hydrophones' are designed so that 
the sensing part of the hydrophone oscillates in the medium 
under the influence of the sound pressure gradient. That 
1s, a hydrophone is constructed for which the driving force 
and the hydrophone’s output is proportional to the 
difference between the pressures at two different locations 
in an acoustic sound field. Therefore, the electrical 
output is proportional to the pressure gradient in the sound 
field. The main characteristic and advantage of a pressure 
gradient hydrophone is the directivity pattern obtained from 
such a small size hydrophone. As will be shown below, 
pressure gradient hydrophones have dipole, or figure-eight, 
directivity patterns, as indicated in Figure 2.13; hence they 


are bidirectional. Assuming the size of the hydrophone is 
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small compared to the acoustic wavelength rr; the dipole 
response when oriented at any angle 4 relative to an 
incoming plane pressure wave is proportional to cos). 

Both for simplicity and since they resemble the fiber 
optic system under study, pressure gradient hydrophones 
formed by combining two small, closely spaced, individual, 
pressure hydrophones will be discussed. The performance of 
this hydrophone pair (or doublet) is easily calculated, and 
the behavior thus derived applies also to pressure gradient 
hydrophones of more complicated construction as long as the 
assumptions of compactness made in the derivation are not 
violated. 

As an introduction to the theory of this type of 
hydrophone, consider two small pressure hydrophones placed 
in a standing acoustic wave field, P(x,t), separated by a 
distance d ced as shown in Figure 2.2. The hydrophones are 
placed so that a null in the standing wave pattern occurs 
halfway between the two hydrophones. The electrical outputs 
of the two are connected 180 degrees out-of-phase with each 
other to obtain ae dipole. The dimensions of the two 
hydrophones are typically much less than a wavelength, 
therefore the acoustic field is assumed to be negligibly 


disturbed. 
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Consider the plane sinusoidal wave shown in Figure 2.3. 


The instantaneous acoustic pressure, P(x,t) is given by: 


i@t 
P(x,t) = P sinlkxle 
o 


(2.1) 


where P is the peak acoustic pressure, k is the propagation 
wave Fees k = Te and x is the distance of one of the 
hydrophones from absolute minimum pressure (i.e. the 
standing wave pressure node). Assuming the argument of the 
sine function is small, the first order (linear) 
approximation (i.e., sin y = y) of the pressure amplitude at 
x is given by: 
jiWt 
Pix,t) =P kx e 


Oo 


(2.2) 


If one sensor is at x = -d/2 and =the other at x = +d/2 as 
shown in Figure 2.3, then the amplitude of the pressure 


difference to be measured may be expressed as: 


Ap =p - P 


: +d/2 -d/2 


(2.3) 
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This may be written as: 


Ap = 


© 


Or: 


Ap = 


For the more general one dimensional case of 


along a standing wave, 


P ktd/Zz — 


(-d/2) J] 


(2.4) 


P kd 


(2.9) 


any point 


consider a standing wave in a tube of 


water as shown in Figure 2.4. The instantaneous acoustic 
pressure P(x,t) is given bys 
ji@t 
P(x,t) = P sin (kxJe 
oO 
(2.6) 
where x has been chosen to be zero at the air-water 


boundary. 
If we expand equation (2.8) 


xX =h, where h is the distance 


the following equation results: 


P(x,t) = P(h,t) 4 OPth,tN\x - 


O x 
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into a Taylor’s series about 


from the air-water boundary, 


WF hc: n+... 
Gf 


h) 


2 


(2.7) 





pressure 
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Figure 2.4 Pressure and Velocity Distribution ina 


Standing Wave Tube. 
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where the first term and all their partial derivatives are 
evaluated at x = h. Time t is held constant. 


The first partial derivative (pressure gradient) iss: 


jWe 
OpsOx = P kcos(kx)e 
oO 
(2.8) 
The second partial derivative is: 
2 2 2 jWet 
C P/Ox = -P k sin(kx)de 
oO 
(2.9) 


The pressure difference between two points axially 
displaced from one another by distance A\x where [Ne =x -h 


is: 
A\p = Pix,t) - Pthytd 
(2.10) 
or 
Ae = «desdxAx + (Oe ™“ rae eee 
(2.11) 


If iv is small compared to the wavelength of the standing 
wave all terms above second order may be neglected and 
therefore the series may be truncated after the second order 


term. 


2 


If h and Ax are selected so that a pressure node is 
located at Ax/2, then the difference in pressure Ape 
reduces to equation (2.5) with Ax =z=d, in the linear 
approximation. The first correction, € , to the linear 


approximation, in the vicinity of a pressure node, becomes: 
€ = -c¢k/\x) /2]tan (kh) 
(2.12) 


A maximum in the pressure difference will occur’ when 
each sensor is equally distant from the pressure node and 
conversely a minimum in the pressure difference will occur 
when each sensor is equally distant from a pressure maximum 
(antinode). 

A pressure gradient hydrophone is sometimes called a 
velocity hydrophone, the implication being that it responds 
to the particle velocity in the acoustic medium. This is 
somewhat misleading since the velocity of the hydrophone 
surface is normally much smaller than the particle velocity 
of the incident acoustic wave. But even though the 
hydrophone is not a direct detector of particle velocity, 
its response, being proportional to pressure gradient, is 
also proportional to particle velocity of the incident 
acoustic wave as shown in the one dimensional = § linearized 


Euler equation: 
-Op/Ox = iWPu (2.13) 
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where @ is the angular frequency of the acoustic wave, () 
the mass density of the medium, andu is the Eulerian 


particle velocity of the free-field acoustic signal. 
B. CALIBRATION OF GRADIENT HYDROPHONES 


Pressure gradient hydrophones may be, and usually are, 
calibrated in terms of pressure rather than velocity 
(pressure gradient is never used as areference). Since the 
relationship between pressure and velocity depends on the 
wave and boundary conditions, these conditions must be 
known, be standard, or fit the definition of sensitivity. 
Therefore the sensitivity of pressure gradient hydrophones 
usually is given in terms of volts/micropascals (V/ [{ Pa) 
specified at a particular frequency. 

Plane progressive waves are specified in the definition 
of free-field voltage sensitivity, and the relationship 
between the pressure p and velocity u in such waves is given 
by the wave or characteristic impedance p/u = Pc- Thus, the 
sensitivity of a pressure gradient hydrophone expressed in 
terms of pressure differs from the sensitivity expressed in 
terms of velocity by the constant factor Pcs where p is the 
density of the medium and c is the speed of sound. 

The standard procedure for pressure gradient hydrophone 
calibration is by direct comparison with a standard pressure 


hydrophone only if essentially plane waves impinge on both 
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hydrophones. Since waves in a free field never are 
perfectly plane or if a spherically diverging wave from a 
point source is used, a correction factor , 
10logli+ eee where r is the distance from the 
source, must be subtracted from the measured sensitivity 
CRef. 71. This correction factor decreases as r increases. 
Because of the difficulties in obtaining free-field 
conditions particularly at low frequencies, a standing-wave 
tube like that shown in Figure 2.4 was employed for the 
laboratory investigations described in this report. Since a 
water-air boundary is a near zero impedance boundary, 
essentially complete reflection will take place at the top 
of the tube. A standing-wave system then can be established 


and the following relations hold in the ideal case (i.e. 


SWR= OO ) 3: 


p =p sin kh, 


oO 
(2.14) 
u= (p (Dc) cos kh, 
oO 
(2.15) 
‘ p/u =Pc tan kh. 
(2.16) 
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If the pressure is measured at a point a distance h from 
the water-air surface, equations (2.14) and (2.15) can be 
used to compute the pressure and velocity at any other 
point. As stated earlier the hydrophones are assumed to 
have a negligible effect on the standing wave pattern. If 
the tube is operated below cutoff for non-axial modes, then 
the pressure variations normal to the tube axis can be 
neglected. A detailed discussion of the calibrator used is 
given in Chapter III section D7. 

A more detailed explanation of sensitivity and 
calibration procedures can be found in Bobber (CRef 7]. In 
this experiment sensitivities of single coil fiber optic 
acoustic hydrophones are defined in terms of HM radians//[/ Pa. 
Sensitivities of gradient (dual coil) hydrophones are 
defined in terms of optical phase shifts induced by acoustic 
Pressures and are expressed in ([lradians/[[Pa/cm. The 
procedures to obtain fiber optic hydrophone sensitivities 
(single and dual coil designs) are discussed in Chapter IV. 
A thorough presentation on fiber optic acoustic sensitivity 


is given by Davis, et al., CRef. 4]. 


C. FIBER OPTIC ACOUSTIC SENSOR 


1. Acoustically Induced Phase Modulation 
Laser light conducted by optical fibers submerged in 
a liquid medium may be modulated (intensity or phase) by 


acoustic pressures. This phenomenon is exploited to make an 
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acoustooptic sensor of underwater sound by demodulating the 
light using a specific detection scheme. Carome and 
Satyshur CRef. 8] present in detail general fiber mode 
theory and the theory of phase and intensity modulation 
induced in optical fibers. Phase modulation occurs when the 
instantaneous phase angle of the lightwave in the optical 
fiber is varied proportionally in accordance with the 
instantaneous value of a modulating signal, e.g. a 
sinusoidally varying acoustic field. 

Two major detection schemes for interferometric type 
sensors discussed in the references are homodyning and 
heterodyning CRef. 5S and 6]. Homodyne detection in optical 
interferometry implies that the reference beam has the same 
optical frequency as the sensing beam. On the other hand, 
heterodyne detection is based on the use or mixing of two or 
more frequencies, i.e@., the frequency of the reference and 
sensing beams differ by LNGe Interferometry is used in 
either case to recover the acoustic signal frequency. 
Homodyne detection was used throughout this experiment so 
only it will be considered in detail. 

An external pressure field applied to a fiber will 
cause various changes in the fibers physical 
characteristics, such as changes of the core radius and 
length, and changes in the values of the optical indices of 
refraction of the core and cladding CRef. 3 and 4]. For the 


case of pressure induced changes of index and of length, the 
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Phase shift A, produced by a change of pressure Ap might be 


given as! 


A=k PRCCOn/GP) + n/XcOPsdr3 
oO 
(2.17) 


where nm is the optical index of refraction of the core, a 
is the propagation constant of light in the fiber, P is the 
acoustic pressure, and 0 is the length of fiber subjected to 
the pressure. The partial derivative (Ons OP) is evaluated 
at constant entropy. Whatever the form of the bracketed 
term of equation (2.17), the phase shift amplitude A enters 
the analysis in identical ways so that it is possible to 
proceed without being concerned with its specific form. 

Using a single frequency coherent laser source, the 
time variation of the electric field vector of the 
lightwave, propagating in a fiber that has been subjected to 
an acoustic pressure field, may be expressed as: 

— = 

E(t) E explit(W t + Asin(W t)]I 

oO oO = 

(2.18) 

where G) is the angular frequency of the coherent laser 
source, Ww nacre angular frequency of the sound field, and A 
is the a tiatis shift amplitude as described above. 


Phase modulation cannot be directly detected by 


Photodetectors due to the fact that they cannot follow the 


S1 


instantaneous values of such high frequency variations 
(approximately vam Hz). Thus, in order to accomplish phase 
detection, an interferometric technique must be used to 
convert the phase modulation to intensity (amplitude) 
modulation prior to detection. A monochromatic laser beam 
is split, with part of the beam being transmitted into a 
reference fiber. The other beam is sent into a sensing 
fiber exposed to an acoustic field. Pressures exerted on 
the sensor fiber as aresult of the sound energy density 
produce optical path length changes in the fiber which, in 
turn, cause phase modulations in the lightwave. After 


passing through the fibers, the two beams are then 


recombined. Thus the total electric field is given by: 


as a <a 
FE 2 E& (t) +E (t) 


T 1 2 
(2.19) 
where 
—_h —* 
E (t) = E explji(G@ t + Asin(W t)] 
1 1 oO s 
(2.20) 
and 
pore § a 
E(t) = E expli(Wt+D) 
2 2 oO 
(2.21) 
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where D is an arbitrary phase difference which will vary 
with temperature, the optical path difference, and other 


— 
factors. E (t) is the electric field vector of the sensing 


1 
a= 
arm and E (t) is the electric field vector of the reference 
2 
arm. The intensity (I) of the recombined beams is 


see 
proportional to the magnitude of the square of E. 
+t 
Therefore, after mathematical manipulation as shown by 


Carome and Satyshur CRef. 8] and neglecting terms that vary 
at angular frequency () and 2G), since as stated earlier 
oO 


o 
they are undetectable by a photodetector, the following 


equation results: 


2 2 amigos 
iy ee ey, Hee Eco J (A) 
gf +5 + EE coms, 
—. ama 
+ 2E- E,sindy (A)sinG)t + 26 E, cosJ, (a)cos2U)t 
+ 2E-E ds (A)sin3WJt + 
e sin sin oe @ 
2 3 Ss 
(2.22) 
— > od 
where E*E makes explicit the fact that the polarization 
tL 2 
directions may not be the same. The first three terms are 


considered the DC terms and all other terms of the series 
are the AC terms. Thus the resulting intensity function 
consists of a series of harmonics of the acoustic 
frequencies where the amplitude of each successive harmonic 


is a function of the acoustic pressure, and varies as the 


Bessel function of corresponding order. 
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If both of the beams are phase modulated by the 
acoustic pressure as in the case of the pressure gradient 
hydrophone, with phase shifts of A and A , respectively, 
the final expression for the otis = identical to 
equation (2.22), except that for A, the quantity 
A\a = A- A must be substituted. 

idee le lightwave with its varying intensity 
is then projected onto the surface of a photodetector. 
There results in the photodetector current a signal which is 
related to the acoustically induced phase shift. The 
photodetector current has components of the following form: 
i(t) =i cos(* {J (k x) + 2) 3 (k ee on 

Oo Oo oOo n=] 2n oOo 


=~i jninp* {2 ye (k x) sin (2n+1) G)t)]I}, 
2n+1 oa 5 


(2.23) 


where J is the Bessel function of order n [CRef. 81, is 
n 
the static Phase difference, which will vary with 


temperature, path length difference, and other factors, 


k =2 Tian >» Where k is the optical propagation constant in 

oO o oO 

the fiber, nN is the optical wavelength in the fiber, and x 
o 

is the amplitude of the optical path-length change at @ . 


S 
The J term in equation (2.23) is the DC component 


Oo 
associated with the acoustically induced signal with the 
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remaining terms being the AC component. Only the AC terms 
of equation (2.23) will be considered in succeeding 
sections. 
2. Feedback Stabilization 

The purpose of homodyne detection in the fiber optic 
interferometer discussed above is to transform the optical 
output of the interferometer into an electrical signal 
proportional to the amplitude of the relative phase shift 
produced by a pressure oscillation at frequency W . 
Homodyne signal processing suffers from phase drift due ea 
small temperature and pressure variations. Any random 
change in length in one arm of the fiber interferometer with 
respect to the other arm appears in equation (2.23) as a 
change in the angle D. As PD changes, the amplitudes of the 
AC components of equation (2.235) varies due to the factors 
cos PD and sin®. This unwanted change in the amplitude of 
the detected signal is referred to as signal fading. 
Consider the case in which the outputs of the two arms of 
the interferometer are initially in phase, D is equal to 
zero, and when combined at the output coupler they interfere 
constructively. If now there is a small acoustically 
induced phase modulation e.g. at peak to peak amplitude O.1 
radians retin to Figure 2.5, the photodetector current 
will oscillate about its maximum value at twice the acoustic 


frequency. If the length of one arm increases due to a 


thermal change, the average output intensity of the 
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interferometer decreases, but the acoustically induced 
oscillations will increase in amplitude and approach a pure 
sinusoidal of frequency equal to the acoustic frequency. In 
a practical hydrophone it is essential to maintain the 
interferometer continuously at its maximum sensitivity. The 
sensitivity to phase changes varies as the slope of the 
curve of Figure 2.5. The maximum sensitivity occurs for odd 
multiples of 7[/2 while essentially zero sensitivity occurs 
for even multiples of 7{/2. Therefore to maintain 
sensitivity a 71/2 phase bias is required. This T[/2 bias 
condition is known as quadrature. 

To place the interferometer in quadrature condition 
for maximum sensitivity and in addition provide a means to 
compensate for slow drifts, a phase-locked compensator 
Circuit is used. The circuitry required to provide and 
insure quadrature operation in the presence of low frequency 
drift is shown schematically in Figure 2.6 CRef. 91. The 
Outputs of the two photodiodes are combined in a 
differential amplifier that provides common-mode rejection 
as well as amplification. This is folloawed Dy two stages of 
integration that provide additional amplification. These 
two integrator-amplifiers pass all signals from DC up to the 
highest frequency of interest. The output of the two stage 
integrator-amplifier is applied to a phase shifter located 


in the reference arm of the interferometer. The phase 
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shifter consists of a lead zirconate-—lead titanate 
piezoelectric (PZT) cylinder around which the fiber in the 
reference arm is rather tightly wound. The output of the 
integrator-amplifier is just equal to the voltage required 
to compensate for phase changes associated with the low 
frequency drift and that induced by the acoustic pressure 
signal of interest. The effect then is to produce a phase 
shift in the reference arm almost equal to that in the 
sensing arm, causing the interferometer to remain very close 
to perfect balance, i.@., to phase-lock the system (Ref. 
10]. If the phase were exactly locked there would be no 
output signal from the interferometer. However, there must 
be an error signal or resulting phase modulation at the 
photodetectors in order to have a feedback signal and this 
is effectively 1/G times the modulation that would exist 
without the feedback, where G is the overall gain of the 
feedback circuit. If the system is initially at a bias 
(operating or quiescent) point away from quadrature, the 
effective gain will not be at its maximum value and the 
system will tend to drift toward its maximum gain point and 
therefore toward quadrature. The output signal of the 
feedback system represents the correction signal required to 
keep the shaw difference between the interferometer at 7(/2 
radians. It contains drift noise and signal information. 
The signal out of the compensating circuit is also 


fad through a high-pass filter that has its low frequency 
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limit set at the lowest frequency of interest but above the 
frequency of thermal drift etc. Therefore, the resulting 
output is a band of frequencies corresponding to the 
frequency range of interest. This constitutes the output of 
the interferometric sensor. 

Operational amplifiers ‘(OPAMPS) are used in the 
feedback circuit and complementary metal oxide semiconductor 
(CMOS) components in the reset circuit. The levels of 
voltage that can be applied by these circuits to the phase 
shifter are the order of + 10 volts. On the other hand, the 
range that the piezoelectric phase shifter can accommodate 
is hundreds of volts. Furthermore, in many cases the 
amplitude of the phase drift resulting from temperature or 
pressure changes is much larger than the phase shift that 
would be generated by applying + 10 volts to the phase 
shifter. Thus, the circuitry keeps track of how much 
voltage has been applied to the phase shifter and if the 
limit of the circuit begins to be reached it automatically 
and rapidly resets the circuit back to its mid range from 
which point it restarts the compensation process. This is 
the purpose of the reset circuit indicated in Figure 2.6. 
The phase change associated with a large amplitude slow 
drift is A maaieaae by a number of saw toothed-like small 
amplitude phase changes. Care must be taken in the 
alignment of the compensation circuit to minimize noise, 


introduced during the reset process. 
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III. EXPERIMENTAL APPARATUS 


A. INTERFEROMETER DESIGN AND CONSTRUCTION 


Figure 3.1 shows a schematic diagram of the experimental 
interferometric system utilized in this study. The 
interferometer is modular in construction. As each module 
was completed proper operation was verified. A brief 
description of the overall construction of the 
interferometer is first presented; this is followed by a 
detailed description of each individual part. 

The system is comprised of an all-fiber Mach-Zehnder 
interferometer. A Laser Diode Labs solid state laser was 
a to high numerical aperture single-mode (830 nm) 
Corning fiber. The laser, with the pigtailed fiber, was 
spliced (see section on fiber preparation and splicing) to 
approximately half a meter of ITT single-mode fiber. The 
free end of the ITT fiber was then spliced to an ITT 2 X 2 
all-fiber 3dB coupler. Initial attempts were made to splice 
the pigtailed Corning fiber directly to the coupler fiber 
but with limited success. Differences between melting 


temperatures -of the Corning fiber and coupler fiber resulted 


1 Pigtailing of fiber to laser accomplished by laser 
manufacturer 
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in low strength, high loss splices. For this reason the ITT 
T-1601 fiber had to be used as an intermediate link between 
the Corning fiber and ITT coupler fiber. 

The couplers themselves are all fiber counterparts of 
bulk optical beam splitters and equally split the light at 
the input end of the interferometer. Half the light is 
coupled to one arm =~and the remaining light into the other 
arm. 

Two single-mode optical fibers with approximately equal 
length (estimated path length difference of approximately 2 
cm) form two arms of the interferometer. Both were 
constructed identically; therefore only the construction of 
one will be discussed. Eight meters, from a ten meter 
length of fiber, was wrapped in a single layer onto a 
Piezoelectric (PZT) cylinder. At each end of the wrapping 
approximately a meter of lead was left. One lead was then 
spliced to fifteen meters of fiber of which approximately 
ten meters was coiled into a 4 cm diameter coil. The other 
lead coming from the fiber wrapped around the PZT is’ fed 
through mechanical polarization controllers (Ref. 111]. One 
end is then spliced to one of the output leads of the input 
coupler and the other end is’) spliced to one of the input 
leads of ene: output coupler to complete construction of an 
arm of the interferometer. The output coupler recombines 
the light from the two arms. Basically the interferometer 


is used to convert phase modulation to intensity modulation 
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which occurs when light from the two arms is combined. The 
two outputs of the interferometer are detected by two 
photodiodes. 

Figure Sed is a photograph of the completed 
interferometer. The design of the system allows for a 
single coil to be placed in water to act as a pressure 
hydrophone or for both coils to be placed in water for 


operation as a pressure gradient hydrophone. 
B. INTERFEROMETER COMPONENTS 


1. Laser Source 

The optical source used in this experimental 
interferometric system was a single frequency infrared laser 
diode optimized to emit coherent light at a wavelength of 
830 nanometers. Two solid state laser diodes from Laser 
Diode Labs were used. The first was a CW-234 laser, 
Pigtailed by the manufacturer with a short section of high 
numerical aperture single mode (830 nm) Corning fiber. The 
second laser, a CWX-233 laser also, pigtailed with Corning 
fiber, replaced the first laser after an inadvertent 
increase in current destroyed the CW-234 laser. Power for 
the laser was provided by a 12 volt 18 amp-hr Gel cell 
battery fhccude a voltage regulating circuit. This circuit 
maintained constant current to the laser diode. Figure 3.3 
is a schematic of the power supply. Before using either 


laser in the interferometer the end of the Corning fiber was 


a 





mi Curews . 2 Photograph of Experimental 


Interferometric System 
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cleaved and set up to enable a Coherent 212 optical power 
meter to sense the laser diode power output. As current to 
the laser diode was increased from O to approximately 53 
milliamps for the CW-234 laser, and O to 39 milliamps for 
the CWX-233 laser, the power of light emitted from the laser 
and coupled into the fiber was read on the power meter. 
Figures 3.4 and 3.5 are graphs of laser power vs source 
current for the CW-234 and CWX-2335 lasers, respectively. 
Extrapolating the linear region back to zero power yields a 
threshold current of 48.5 milliamps for the CW-234 laser and 
31.3 milliamps for the CWX-233 laser. The current to the 
laser diode was maintained at 53 milliamps for the CW-234 
laser and 33 milliamps for the CWX-233 laser during 
experimental testing and data acquisition. 
2. Fiber Specification 


The optical fiber used was ITT single-mode fiber, 
Ps 
type 17-1601 optimized at 0.83 microns. Specifications of 


fiber are as follows: 


Fiber Ident: 81092U-17Bic 
Preform No. : EMT-21972B 
Numerical Apertures 0.12 

Goce Disseter: 4.57 microns 


2 Optical fiber obtained from International Telephone 
and Telegraph, Electro Optics Products Division, Roanoake, 
Va 24019 
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Outer Cladding Diameter: 75 microns 


Primary Sheaths GE 615 Silicone 

Secondary Sheaths polyester Hytrel 

Total Diameter: 16 mil 

Attenuation: 2-907 dB/Km at 0.83 microns 


Figure 3.6 is a sketch, provided by the manufacturer of the 
end of the fiber, with typical dimensions, approximate index 
of refraction profile, and typical spectral attenuations. 
S- 3 dB Coupler 
The 2X 2 single-mode couplers were also purchased 


from ITT. The following specifications are provided: 


Serial Nos.: JM-SM-108 & JM-SM-109 
Fiber No.: EMC 41556C/830427-401a 
Fabrication Date: 6/3/83 

Excess Loss: O.1 dB 

Uniformity: O.1 dB 

Operating Wavelengths: 0.83 microns 


The purpose of a3 dB coupler is to split the light equally 
into the two arms of the interferometer or to recombine the 
light causing the light from the two arms to interfere in 


the output fibers and on the face of a photodiode. 
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4. Piezoelectric Phase Shifter 


The phase shifter consists of a lead zirconate-lead 
titanate (PZT) eat around which the fiber is tightly 
wrapped. The cylinders used were type Channel S500. Each 
cylinder is 3.8 cm (1.5 in.) long by 3.8 cm in outer 
diameter with a wall thickness of 0.32 cm (0.225 in.). Two 
of these cylinders were incorporated into the 
interferometer, one in each arm. The PZT cylinders have a 
relatively small dynamic range but by wrapping many turns of 
fiber on the PZT it was possible to increase the magnitude 
of the phase shift. Chapter II, Section C presents’ the 
theory of operation of the PZT as a phase shifter. 

3- Polarization Controller 

A polarization controller, as described by Lefevre 
CRef. 111, was used in this experiment. A very simple 
configuration is employed to control the polarization of 
single mode fiber. This device is equivalent to fractional 
wave plates of classical optics. The controller uses the 
stress birefringence induced by bending the fiber. 

6. Photodetector 


The photodetector used to detect the optical output 
4 
of the fiber interferometer is a Clairex type CLD-41 


3S Piezoelectric cylinders obtained from Channel 
Industries, Inc., Santa Barbara, California 93105 


4 photodiodes obtained from Clairex Electronics, Mount 
Vernon, New York, New York 10350 


w2 


photodiode. They are all silicon PN planar diodes with high 
linearity, low dark current, and fast response. The 


following are the electrical characteristics as given by the 


manufactures 
Active Areas 1.3x1.3 mm 
Short Circuit Current: min 6 to max 12 microamps 
Open Circuit Voltage: -40 volts typical 
Dark Currents 1 microamp 
Junction Capacitance: 200 picofarads 
Rise or Fall Time: 2 microseconds 


Temperature Coefficients +.24/degree C typical 

Peak Spectral Response: 0.91 microns 
Figure 3.7 is a plot of photodiode DC voltage vs’ system 
laser optical power as measured in the laboratory using a 
Hewlett-Packard HP-S438A digital multimeter to measure the 


output voltage. 
C. FIBER PREPARATION AND SPLICING 


Before splicing between fibers could be accomplished, 
each end of the fibers had to be prepared by stripping and 
Cleaving. Fiber stripping was accomplished in two steps. 
First, the outer Hytrel plastic sheath was removed by use of 
a razor blade. The fiber was placed on a piece of black 
construction paper and then a razor placed at a very small 
angle to it. With a slight pressure applied to the razor, 


the fiber was drawn out from underneath the razor. This 


re 
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will carve the Hytrel jacket neatly from the fiber. By 
rotating the fiber and repeating this process several times, 
most of the outer Hytrel and inner silicone sheathing 
material will be removed. 

The second step, that of removing any residue of inner 
or outer sheath that was not removed by the razor blade, is 
best done chemically with a solution of chromic acid. This 
is done by immersing the section of fiber that has been 
partially stripped by the razor blade in the acid (for 
approximately two minutes. Any jacket material that remains 
will turn to a jelly-like consistency which is’7 easily 
removed by wiping with a soft tissue. After rinsing in 
distilled water, the fiber is pulled between a _ tissue 
saturated with methyl alcohol. If the fiber is completely 
Clean of the jacket a squeaking sound will be heard. At 
this point the fiber is bared to the silica cladding and 
must be handled with caution. An unscratched fiber has very 
high tensile strength, but will break easily if scratched 
and bent. 

In order to insert and extract light efficiently to 
obtain a low loss splice, the ends must be flat = and 
perpendicular to the longitudinal axis of the fiber. The 
process utilized to achieve this was cleaving. To cleave 
the fiber, the ends are first stripped back about two 
centimeters longer than needed. After preparation for 


Cleaving, the fiber was held with slight tension and 


raf 


carefully touched to the Cleaving blade at a point 
approximately 2 cm from the end of the Bain The cleaved 
area was then placed on a smooth, round mandrel that has a 
radius of approximately four centimeters. The fiber was 
held with thumb of one hand. Maintaining modest tension, 
the fiber was pulled around mandrel with other hand. When 
scored area reaches mandrel, acliean break will usually 
result. Experience showed that if the fiber broke while 
being scored it would have a bad end face. 

To check the cleave, the fiber is inspected under a 
microscope. Angular breaks other than 90 degrees will not 
make a good low loss splice. In addition, chips and barbs 
at the cleaved ends will not allow satisfactory splices. 
Figure 3.8 shows a good cleave along with improper cleaves. 

Splicing of fiber was accomplished using an electrical 
arc welder. The type of welder used was a Power Technology, 
Inc. portable fiber splicer model PFS-240. Basic splicer 
operations of prefusing and fusing for monomode fiber were 
done in accordance with the operating manual. Thumbwheel 
settings which controlled splicer arc amperage and the time 
variation of the arc current usually used are as follows: 

Ramp Time Arc Time Arc Current 


- 


0.1 sec 0.6 sec 16 mA 


S Sapphire cleaver obtained from General Fiber Optics, 
Caldwell, New Jersey 07006 
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Final alignment of cores for maximum light transmission 
was used to produce a low-loss splice. The 
micro-positioners on the splicer were adjusted to bring the 
prepared ends of the fiber into close alignment. Light from 
the laser was injected into the free end of one of the 
fibers, transferred across the proposed splice to the second 
fiber, and the output from the latter was projected onto the 
face of the power meter’s photodiode. The micropositioners 
were adjusted until the power meter indicated that the 
transmission across the splice was maximum. At that point 
the fiber ends were brought together, heated by the arc 
welder and fused. As the joint cools, surface tension will 
pull the fibers into final alignment. Reports in the 
literature indicate that losses obtained by fusion splicing 
may be as low as 0.1'’s dB. However, in this study many 
hours of practice were required to obtain good (i.e. less 
than 3dB loss) splices. Nine splices are required in the 
configuration of the experimental system described above. 
The average loss from each splice was approximately 1dB and 
the total loss of the constructed interferometer was 
approximately 8.4dB. 

Prior to projecting light from the end of the fiber onto 
the irraree Of a photodiode, the end of the fiber was 
prepared in the same manner as for splicing with the 
exception that the length of fiber that had its sheath 


removed was painted black. Painting was used to accomplish 
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mode stripping. Cladding light was stripped prior to 
reaching the end of the fiber, thus only core light was 
detected by the photodiode. In splicing each additional 
fiber element into the system it was necessary to mode strip 
its output end just before the energy meter to ensure that 
only core light is detected, i.@. to absorb unwanted light 


that may reach the output end via the cladding. 
D. INSTRUMENTATION AND DATA ACQUISITION 


Basic block diagrams of the experimental systems used to 
obtain data on interferometer performance will be shown in 
Chapter IV as the experimental procedures are discussed. 
Computer data acquisition was used for portions of the data 
taking. The Hewlett-Packard HP-85F computer coordinates the 
peripherals and records and displays the data. The printer 
and plotter are non-essential to the system and were used to 
provide larger format graphics and printouts than are 
available with the built-in thermal printer provided with 
the HP-8SF. The following is a brief description of each 
individual piece of equipment. 

1. HP-85F Computer 

The HP-SSF is an eight bit microprocessor that 
utilizes BASIC computer language. The computer has as 
standard 16K bytes of read/write memory. 16K bytes of 
additional memory were added to the system giving a total of 


S2K bytes. The HP-SSF has a 127 millimeter diagonal black 
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and white electromagnetic CRT. A 32 character per line 
thermal printer/plotter is part of the unit. Programs or 
data may be stored on and read _ from magnetic tape 
cartridges. To interface with peripheral equipment, an I/0 
ROM and an interface card were added to provide HP-IB (IEEE 
standard 488-1975) instrumentation capabilities. 
2. Arbitrary Waveform Generator HP 33514A 

The Hewlett-Packard Model 3314A arbitrary waveform 
generator can produce three different pre-programmed 
waveforms} sine, square, and triangular. This instrument 
has a frequency range of 1 =microHertz to 20 megaHertz. 
Output amplitude is O millivolts to 10 volts peak to peak 
into a fifty ohm load. This model is fully programable 
through a HP-IB connection. For this experiment sine waves 
of various frequencies and amplitudes were usually 
generated. During alignment of the locking circuit, 
however, a small amplitude 100 Hertz square wave was used. 

S.- Spectrums Analyzer HP SS82A 

The Hewlett-Packard Model S3582A is a dual channel 
spectrum analyzer. This instrument has a frequency range of 
0.02 Hertz to 25,600 Hertz. The analyzer has a 11.9 by 9.6 
cm CRT that can display two simultaneous information traces, 
plus four lines of alphanumeric data giving measurement 
configuration and results. Frequency spans from 1 Hertz to 
29,9000 Hertz full scale allow flexibility in selecting the 


portion of the spectrum to be analyzed. Spans’ from 5 Hertz 
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to 25,000 Hertz can be positioned anywhere within the 
frequency range of the instrument to provide excellent 
frequency resolution. The instrument's "front-end" 
sensitivity ranges can measure and analyze from 1 micro-volt 
to 31.6 volts and has a dynamic range of 70 GB. 
4. Oscilloscope Tektronix 2215 
The Tektronix Model 2215 oscilloscope is a dual 
channel, 60 MHz instrument. Its vertical system provides 
calibrated deflection factors from 2 millivolts to 10 volts 
per division. Trigger circuits enable stable triggering 
over the full bandwidth of the vertical system. The 
horizontal system provides calibrated sweep speeds from 0.5 
seconds per division to 50 nanoseconds per division. 
3. Digital Multimeter HP 3438A 
This multimeter was utilized to obtain rms voltage 
readings of an LC-10 hydrophone output and transmit this via 


the HP-IB to the HP-85 computer. 


6& Standard Hydrophone LC-10 


6 
An LC-10 hydrophone , Serial No. 2341 was used as 


the standard hydrophone for sensitivity calculations of the 
single coils and of the dual (gradient) coil setup. Average 
free field voltage sensitivity for this hydrophone is 


specified by the manufacture to be -209.2 dB re 1 Volt/[ Pa. 


& Hydrophone manufacture is Celesco Ind., Canoga Park, 
California, 91304 
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This hydrophone and an USRD type F27 reversible transducer 
were used in another experiment (Ref. 12] to obtain a 
reciprocity calibration of an USRD type F33 transducer. The 
sensitivity for the LC-10 hydrophone obtained from the 
experiment was within 2 dB of manufacturer's specifications 
between 5 an 3O KHz. 
7. Acoustic Calibrator 

An acoustic system was constructed for calibrating 
hydrophones of various types by the comparison method. 
Figure 3.9 shows the calibrator in the lab, and Figure 3.10 
shows the details of construction. The calibration section 
consists of an open-ended column of water excited by a 
moving-coil driver. A USRD type Jil transducer is used as 
the driver. The calibration chamber consists of a & inch 
inner diameter, 1/4 inch thick PVC 1120 tube (ASTM 1785), 
mounted on aeorim around the face of the driver. 
Compensation for the hydrostatic pressure of the water on 
the face of the driven piston is provided by the airtight 
enclosure mounted under the driver assembly. A manometer is 
used to measure the equalizing pressure. A valve is opened 
and air is pumped into the equalizing chamber until the 
internal air pressure is equal to the water pressure on the 
upper face of the driver. This system is similar in design 


to the USRD type G40 shipboard calibrator used by the U.S. 


Navy CRef. 131, based on article by Schloss and Strasberg 
CRef. 14]. 
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Figure 3.10 Sectional View of Acoustic Calibrator. 


The compliance of the PVC tube reduced the speed of 
sound of the water in the calibrator to approximately 1/5 of 
the normal speed of sound in water. This was useful in the 
experiments since it reduced the wavelength IN by a factor 
of S at any given frequency and resulted in a larger 
Gradient. Figures 3.11, 3.12, and 3.135 show the standing 
wave acoustic field at three resonant frequencies of the 
calibrator. The figures are graphical plots of LC-10 
hydrophone voltage as a function of position in the 
calibrator. By measuring the distance of a half wavelength 
at the three resonant frequencies the speed of sound in the 
calibrator can be calculated. From the data shown in Table 
I an average value of 358.3 meters per second with a 
standard deviation of 8.6 m/sec for the speed of sound in 
the calibrator is calculated. The data was gathered over a 


period of one month. 


TABLE I 


Sound Speed in Calibrator 


Frequency (Hz) Hal f Speed of 
Wavelength (cm) Sound (m/sec) 
390 44.5 347.1 
390 44 343.2 
391 435.5 355.8 
375 390 357.0 
39S . 1 367.7 
39S 30.5 361.7 
IIS 30.35 361.7 
39S 30.35 361.7 
7352 24.5 368.5 


average >>> 358.3 + 8.6 
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8. Locking Circuit 
A phase lock loop locking circuit similar to the 
circuit presented by Fritsch and Adamovsky [CRef. 9] and 
discussed in Chapter II was used to verify operation of the 


interferometer as a single and dual (gradient) hydrophone 


system. Data will be presented in the next chapter. 
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IV. EXPERIMENTAL PROCEDURE AND RESULTS 


A. INTERFEROMETER CHARACTERISTICS 


After construction of the interferometer was completed, 
testing was conducted to verify proper operation of the 
interferometer as a sensor. The initial tests consisted of 
detailed measurements taken to determine the amplitude of 
the optical phase shift as a function of the drive voltage 
applied to the piezoelectric cylinder in one of the arms of 
the interferometer. Figure 4.1 shows a block diagram of the 
instrumentation required for operating the system. A 
Variable amplitude sine wave of frequency f from the 
function generator was applied to the PZT cylinder. & 
Tektronix 2215 oscilloscope and a HP 3582A spectrum analyzer 
were used to detect the output of one of the photodiodes. 

The behavior of a Mach-Zehnder fiber interferometer was 
treated in detail in Chapter II. Recall that the AC portion 
of the Photodiode current is given by the following 


equations 


i(t) = i soe * (k Ncostan( wey 


OO 
2D 3, 
—i sti [eos Ck Wine (2ned) (GW) £93] 


O 2n+1 oOo Ss (4.1) 


where J is the Bessel function of order n. 
n 
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Figure 4.1 Block Diagram of Instrumentation Setup 
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From equation (4.1) it is clear that as the phase shift 
k x is increased the photodiode output becomes non-linear 
a harmonics of the drive signal begin to appear. This is 
illustrated in Figure 4.2 where oscillograms of the 
photodiode output signal are shown for successively higher 
drive signals applied to the piezoelectric cylinder. In 
each oscillogram the drive signal is shown in the upper 
trace and the corresponding photodiode output signal is 
shown in the lower one. 

The magnitude of the optical phase shift as ae function 
of drive voltage may be determined quantitatively by 
examining the frequency spectrum of the photodiode output 
signal. From equation (4.1) it may be seen that the 
amplitudes of the various a order harmonics should have 
nulls for phase shifts kx equal to the zeroes of the 

th ° 
corresponding n order Bessel functions. This behavior is 
illustrated in Figure 4.3 where the photodiode voltage is 
plotted for the fundamental and first two harmonics vs drive 
voltage applied to the piezoelectric cylinder. The data 
were obtained with the computer data acquisition equipment 
discussed in chapter III. A computer program for the HP-8S5F 
computer, listed in Appendix A, was written in BASIC 
computer language to control the peripheral equipment. A 
fundamental frequency, lower and upper limits of amplitude, 


and an incremental value are selected. The function 


generator and spectrum analyzer are preset and configured to 
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the prescribed settings. At the lower voltage amplitude 
limit the spectrum analyzer seeks for the maximum amplitude 
of the fundamental and first two harmonics attained during a 
selected period of time. The amplitude of these signals 
vary in time as the angle PD varies randomly. The sample 
time selected on the spectrum analyzer, of the order of 14 
seconds, was sufficient that the actual maximum values were 
attained and recorded. The HP-835F computer directs the 
spectrum analyzer to output the maximum amplitude values of 
the fundamental and harmonics. These values are stored by 
the computer. The computer then directs the function 
generator to increase the drive voltage by the specified 
incremental value. The spectrum analyzer again seeks for 
the maximum amplitudes and the process is repeated. This is 
done for each incremental value until the upper limit is 
reached. The data stored in the computer memory is sent to 
the printer and also transferred to magnetic tape for 
permanent storage to enable analysis at a later time. a 
Graph of the data is plotted on the computer. printer. 
Sample of data obtained is seen in Appendix B. From this 
data plots of the type shown in Figure 4.3 could be 
gener ated. 

The computer program was run initially with the CW-234 
laser source. Three runs were made, two with the fiber 
wrapped PZT cylinder of one arm acting as the modulator and 


one with the other fiber wrapped PZT of the other arm. 
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After replacement of the source with the CWX-233 laser was 
accomplished two more runs were made, one on each PZT 
cylinder. From the data, the amplitude of the drive signal 
that causes the fundamental to be at its first null was 
determined to be 1.65V (peak) +0.04V(peak). Thus, referring to 
equation (4.1), it may be concluded that the drive condition 
was such that the first order Bessel function J was ata 
zero. The first of the fundamental’s zeroes Baloo! at a 
phase shift of k x = 3.83 radians. Therefore for the 8 
meters of fiber faa on the PZT cylinder a value of 0.29 
rad/volt (peak) /meter results. By observing the drive 
voltages required to produce the first zero of the 
fundamental as the frequency was varied manually it was 
possible to determine that the phase shift per volt remained 
constant, within the precision of the measurements, for 
frequencies from 80 hertz to 1500 hertz. 

Agreement of experimental data presented in Appendix 8B 
and illustrated in Figure 4.3 to Bessel function theory is 
shown in Table II. The ratio of experimental maxima and 
minima of the same harmonic are compared to the theoretical 
ratio of maxima and zeros of Bessel functions of the same 
order. Theoretical data as tabulated by Abramowitz and 
Stegun [Ref. 15] was used. 

The last test conducted was to verify operation of the 


system with the locking circuit. After aligning the 


locking, a small amplitude sine wave was applied to one of 
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the PZT cylinders simulating an acoustic signal. The 
locking circuit produced a signal, similar to this input 
signal, that was fed back to the PZT cylinder in the other 
arm to act as a compensation signal, thus maintaining the 
two arms in quadrature. This output signal from the locking 
circuit was also sent to the oscilloscope to show that the 


output signal did follow the input signal. 
B. INDIVIDUAL SENSOR SENSITIVITY 


The individual acoustic sensitivity of the sensor coils 
in each arm of the interferometer was calculated from 
measurements obtained in the acoustic calibrator. A single 
coil was mounted on a stand with an LC-10 hydrophone placed 
in the center of the coil as shown in Figure 4.4. The 
acoustic calibrator was driven at a frequency that generated 
a standing wave acoustic resonance. Chapter III section D7 
discusses the standing waves generated in the calibrator. 
The coil/LC-10 hydrophone setup was lowered to the first 
maximum point of the standing wave. Output of one of the 
photodiodes of the interferometer was sent to the Tektronix 
2215 oscilloscope and the HP-3582A spectrum analyzer. 
Output of the LC-10 hydrophone was sent to the HP-3438A 
Multimeter eat to read AC voltage. A schematic of the 
instrumentation setup is shown in Figure 4.5. The spectrum 
analyzer was set to monitor the fundamental frequency and 


the first two harmonics. Drive voltage was increased to the 
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Figure 4.4 Photograph of LC-10 Hydrophone and 
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acoustic calibrator until the voltage of the fundamental 
frequency, as seen on the spectrum analyzer, reached its 
first minimum value. The voltage of the LC-10 was 
recorded. If there was more than one pressure maximum the 
sensor arrangement was moved to the next pressure maximum 
and the LC-10 voltage at the fundamental ’s first zero point 
was recorded. Data were recorded for all of the pressure 
maxima of the standing wave. Data were obtained both 
manually and by computer control as discussed in the above 
section. 

The first minimum of the fundamental occurs at a phase 
shift of 3.83 radians, therefore the sensitivity of the coil 
M in optical radians per micropascal is calculated as 


Cc 
follows: 


M= M #3.835rads/V 

C H H 
(4.2) 
where M is the sensitivity of the LC-10 hydrophone in 
volts/ [Pa obtained from manufacture specification sheets. 
V is the LC-10 hydrophone voltage where the first zero of 
Be fundamental occurred. Table III is a summarization of 
all single coil sensitivity measurements taken at three 


different calibrator resonances and different depths 


associated with pressure maxima. 
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Table III 


Individual Coil Sensitivity 


Coil # 1 Coil # 2 

Freq:Depth of iNos. iAvg iStd 

(Hz) :Hyd from : of | M = idev 
(Air-Water !Runs (Ur ad wir ad 


8 
Freq:iDepth of iNos.iAvg iStd : 
(Hz) i:Hyd from : of | M idev i 

‘Air—-Water ! Runs Uirad (Lr ad 0 


‘Boundary | iL Pai Li Pa} ‘Boundary | i Pai {/ Pa 
i (cm) i bo 93 t 3 3 ; (cm) i ( 3 t -3 
h i (10 110 |¢ ‘ i (10 $10 
392 | 23 ; 2 $7.07: -— $1390 | 23 i 1:15.65! - 
3835 i 60 © 2 $12.4f -— {1389 | 65 ; 1 t9.99i - 
385 i 80 : YF {9.13512.31°°385 | 78 i; 9Y $110.3:12.4 
i i i i tt i i 0 5 
D44% | 16 i 6 $4.95:1.563:':593 | 16 i 4 13.19:.314 
393 | 48 : 7 $8.6611.14:°1593 | 48 ; S$ $11.5'%.529 
344 | 47 f 10 19.93:2.9715 — ¢ = ; - | = tii 
1 i n ‘ ft i i : 5 
695%! 12 : 9 19.77{2.41:11704 | 11 ; 13 1:8.4511.36 
695 | 37 ; § 110.41. 2951:695 | 36 | 18 $12.611.57 
6935 | 65.5 { 7 i125.4:4.991:700 | 62.5 i: 18 121.6:9.15 
6935 } 91.5 { 28 118.616.8 11'690 |} 92.5 [{ 6 1[20.719.27 


* Depth of water column changed from 93cm to 100cm to 
shift second resonance from 393 to 544 Hz, and to 
shift third resonance from 7352 to 4695 Hz. 

If all of the measurements taken are grouped together 


neglecting frequency and coil position in the calibrator, 


the blind average of the sensitivities are: 


3 
Coil #1 M = 13.8 xX 10 Lrad/ [Pa 
Ci 


with a standard deviation of: 


3 
O0 =7.5X 10 [rad/[[Pa 
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and 


3 
Coil #2 M =13.9%X 10 {[rad/[{Pa 
C2 


with a standard deviation of: 
=> 
O=7.7X 10 [Lrad/ [Pa 

As can be seen from Table III and the blind averages, the 
individual coil sensitivities are matched well within 
experimental error. 

Table III indicates an increase in coil sensitivity as 
depth is increased with exception of very deep depths where 
proximity of the sensors to the source appears to have 


effected the measurements. 
C. GRADIENT SENSOR SENSITIVITY 


A gradient sensor was formed by placing the two fiber 
coils in the calibrator separated by a distance of 10 cm. 
The LC-10 hydrophone is) placed in the center of the coils 
approximately 5S cm from each coil. The LC-10 hydrophone and 
both of the coils are then submerged in water and exposed to 
an acoustic field and an established acoustic pressure 
gradient. If the pressure sensitivities of the two coils 
are closely matched, the pair does not respond to the 
pressure, but to the spatial pressure gradient. Such a 


Gradient sensor can be utilized to detect the direction of 
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an incident acoustic wave as discussed under conventional 
gradient hydrophones in chapter II. 

The J11 transducer source of the acoustic calibrator is 
driven at a frequency to produce a well defined standing 
wave. The LC-10 hydrophone/ dual coil setup is submerged to 
a location of a pressure minimum, defined as a point where 
the voltage from the LC-10 hydrophone is at a minimum. With 
the LC-10 at this position the two coils are located on 
either side of the pressure minimum and therefore in a 
region of maximum pressure gradient. The external 
instrumentation arrangement is the same as when § individual 
coil sensitivities were obtained. The driver amplitude is 
increased until the amplitude of the fundamental frequency 
reaches its first zero point. However, the drive voltage, 
not the LC-10 voltage was recorded. Since the LC-10 
hydrophone is at a minimum pressure position, its output 
voltage is quite low and, therefore, the pressure in the 
region of the coils could not be directly computed using 
just LC-10 sensitivity. 

However, by knowing LC-10 voltage (V) at a pressure 
maximum for a set drive voltage a ee between LC-10 
voltage and Jil drive voltage could be calculated. Using 
LC-10 sensitivity (M ), a value of pressure per J11 drive 
volts is obtained. This is shown in the following formula: 


P /Jlilvolts = [€V at P J/[IM *Jlivoltage] 
Oo H max H (4.3) 
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Recall from Chapter II, equation (2.6), the pressure at 
a distance x above and below (+,-) a pressure’ minimum 


(assuming x << d) iss 


P = +P kx 


1+ 
oO 


(4.4) 


Replacing P of equation (4.4) with equation (4.3) and 
oO 


multiplying by the drive voltage needed to reach the first 


zero of the fundamental (V ) the pressure at a distance x 
Jil 
iss 


= C(P /Jlivolts) *#V kx 
o Jil 


140 


(4.5) 


With the coils equally spaced on either side of the 
pressure minimum, the difference in pressure is twice the 
pressure at one of the coil locations. The pressure 
gradient is the pressure difference divided by the 


separation between the coils, and is given by the following 


formulas 


Vp = 2P //Ax 
+ 


(4.6) 
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Gradient sensitivity is calculated directly by the 
following formula, when the pressure gradient is sufficient 
to null the amplitude of the fundamental: 

M = 3.83rads/\/P 
D 


(4.7) 


The units of gradient sensitivity are radians/ Ll Pa/cm. 
Gradient sensitivity can also be calculated indirectly 
by using single coil sensitivity values. Indirect gradient 


sensitivity can be expressed as follows: 
" = (6 -¢ vr 
ID P, Vv 


where D and D are the phase shifts in coil 1 and coil 2 
1 2 


(4.8) 


respectively. Using the pressure value as expressed in 


equation (4.5) the phase shifts may be expressed as: 


P 


ll 

a 
a 
= 


1 + C1 


(4.9) 


and 


(4.10) 
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where P and P are the pressure amplitudes at the sensor 
coil Si icons bat M and M are single coil sensitivities 
obtained from ere aT ie Tl ttes equations (4.9) and 
equation (4.10) into equation (4.8) and using the fact that 
P = -P when the sensor is at a pressure node, the indirect 
ea dot sensitivity can be expressed as: 


M =P CM +M”™ a/\VJP 
ID + Cl C2 


(4.11) 


Substituting equation (4.6) into equation (4.11) results in 
an indirect gradient sensitivity formula as follows: 
mM = tM +m iI/c2//\x3 
ID Ci C2 


(4.12) 


where f\x is the separation between the two coils and this 
distance is assumed to be smaller than a wavelength. 

Table IV is a summary of direct and indirect gradient 
sensitivity measurements calculated from dual coil and 
single coil data respectively. Measurement of direct 
gradient sensitivity was calculated at three different 
calibrator resonances and different depths associated with 
pressure minima (maximum pressure gradient). Indirect 
Gradient sensitivity is also calculated at three different 


calibrator resonances but uses the average sensitivity 
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obtained at a pressure maxima for the single coils divided 
by the desired distance between the coils (equation (4.12)). 

Error analysis for direct gradient sensitivity can be 
made by taking the square root of the sum at squares of the 


constituents of equation (4.6). This may be expressed as 


follows: 
pma a 2 2 
Ee. 
OM\ _ OM) , (Aun) , (0 (OX 
M a a V eo xX 


(4.13) 


where C is the experimental average speed of sound in the 
calibrator as shown in Table I. All other symbols are as 
explained in equation (4.3) and equation (4.5). 


Error analysis for indirect gradient sensitivity may be 


expressed as follows: 
2 2 |/ 
Sit) = | (CF). Ox 
Af x 
der i 


(4.14) 


where <M> is the average deviation for single coil 
sensitivities and M is the average of the single coil 


sensitivities as presented in Table III. 
D. ANALYSIS 


A general analysis of single coil sensitivity compared 
to published data (Ref.6] can be shown. From Table III 


typical sensitivities are on the order of approximately 
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10 [rad/[ Pa per ny GEeP naw ee 10 meters of fiber. This 
results in a sensitivity of 10 /[/rad/[{ Pa/m. For an optical 
wavelength of 0.83 microns, the experimental value of 
A d/pAP is approximately 130 x ie ae al For the 
same diameter fiber (0.4mm), the calculated value is greater 
than published values (Ref. 3] by a factor of four to five. 
This variation could be explained in terms of the large 
variations in observed physical properties of Hytrel 
coatings resulting in a more sensitive fiber. 

As shown in Table III, the single coil sensitivities of 
coil 1 and coil 2 at the same depth appear to be closely 
matched within experimental error. The sensitivity of the 
individual coils are not matched over a depth range in the 
acoustic calibrator. The data indicates that the 
sensitivity of the single coil is a function of depth with 
possible near field effect occuring close to the source 
causing an anomaly. No reason could be found for this 
result but follow-on experiments will investigate this depth 
dependent sensitivity of single coils. 

Measured pressure gradient sensitivity and calculated 
pressure gradient sensitivity as shown in Table IV in the 
previous section are in agreement within experimental 
error. The measured pressure gradient sensitivity obtained 
from data at pressure minima is shown to be more constant 


with frequency and depth than the calculated pressure 


gradient sensitivity based on single coil measurements at 
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pressure maxima. Where the single coil data at Table III is 
a function of depth the measured sensitivity of the dual 
coil setup appears to be relatively independent of depth 
with an average value of 0.124+0.015 [LLrad/L Pa/cm. 

Final analysis of the pressure gradient hydrophone was 
conducted by connecting the output of the photodiodes to the 
phase-lock loop locking circuit. The dual fiber coils with 
a LC-10 hydrophone in the center was lowered in the acoustic 
calibrator. The voltage from the locking circuit and the 
LC-10 hydrophone were recorded. Figures 4.& and 4.7 are 
plots of the recorded data at S93 Hz and 752 Hz 
respectively. As shown in Figures 4.464 and 4.7 as LC€C-10 
hydrophone proceeds to a minimum gradient hydrophone voltage 
proceeds to a maximum and as LC-10 hydrophones proceeds to a 
maximum gradient hydrophone voltage proceed to a minimum. 
The data to produce Figures 4.6 and 4.7 is in good agreement 
with theory except for data recorded near the source or near 


the air-water boundary. 
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Figure 4.6 


Output at 593 Hz. 
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Figure 4.7 


Output at 752 Hz. 
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Ve. CONCLUSIONS AND RECOMMENDATIONS 


Since their invention in the late 1970’s CRef. 3 and 4], 
many applications for fiber optic sensors have been 
suggested. In the area of low frequency (audio) sound 
transduction in the sea (hydrophones), it was felt the most 
reasonable application of this class of sensor was as a 
dipole sensor. The potentially high sensitivity could be 
exploited since the output of a dipole sensor is the small 
difference of two relatively large signals. Secondly, an 
interferometric (homodyne) detector is intrinsically 
differential since, in this case the interference is created 
by the difference of the optical path lengths in the two 
arms of the interferometer. An additional benefit of this 
differential operation is the high degree of common mode 
noise rejection which maybe gained by having the "sensor" 
and “reference" arm in the same environment until the very 
end of the sensor. 

The laboratory investigations outlined in this’ report 
have demonstrated the feasibility of this application. 
Because the fiber optic system is essentially a phase 
modulated eehen than amplitude modulated system (equation 
2-22), it is difficult to make an absolute comparison 
between it and a piezoelectric pressure gradient sensor. 


Before the impedance converting, 20 dB gain, internal 


94 


amplifier, the USRD type G63 standard pressure gradient 
hydrophone CRef. 13] has a sensitivity of 7.5X10 V/[LPa/en. 
The fiber optic system described here using piezoelectric 
phase shifter (Chapter III section A-4) and the feedback 
locking circuit (Chapter III section D-8) had an output of 
305+8 mV/radian (Chapter IV section A and Fig. 4.3). This 
combined with an average gradient sensitivity of 
0.124+0.015 [Lrad/{{Pa/cm TMS overall sensitivity for 
the fiber system of 3.8+0.5X10 V/LL Pa/cm. By this method, 
the fiber system has a 14 dB advantage over the USRD type 
G63. 

At this stage, further work in this area is justified. 
The directions that should be pursued are the ruggedizing of 
the sensor, and the development of a package which will 
incorporate the sensor and couplers ina single submersible 
unit, possibly to be followed by the integration of the 
laser and demodulators in the same unit. Alternative 
demodulation schemes which do not require a piezoelectric 
fiber stretcher to maintain “quadrature lock" should also be 
investigated. Finally, direct comparison with other 


directional dipole sensors such as those on DIFAR sonobuoys 


should be made under real sea conditions. 
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DATA ACQUISITION PROGRAM 


This program was written to use a HP-85SF computer to 
gather data more accurately and quickly than by manual 
means. This program was used after initial construction and 
prior to sensitivity experiments to obtain data for 
computation of interferometer characteristics (see Figure 
4.1). Chapter IV Section A gives a detailed explanation of 
how this program was utilized for computerized data 


acquisition. 
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OPTION BASE 1 

SHORT F(3), A(250), B(250,3), C(250,3), E250) 
REM *#* OBTAIN INITIAL VALUES FOR PARAMETERS ** 
DISP “ENTER FUND FREQ" 

INPUT F(1) 

DISP “ENTER SPECT ANAL START FREQ" 

INPUT Fil 

DISP “ENTER START AMPLITUDE IN MVOLTS" 

INPUT A(1) 

DISP “ENTER END AMPLITUDE IN MVOLTS" 

INPUT V1 

DISP “ENTER INCREMENT VALUE IN MVOLTS" 

INPUT I 

F (2) s2eF (1) 

F (3) =3#F (1) 

REM *#* PRESET AND INITIALIZE HP-3314 AND HP-33582 *+* 
OUTPUT 707 3; "PR" 

OUTPUT 707 3 "“FR",F (1), "HZAP,AC1) , "MVFU1L" 
OUTPUT 711 3;"PRS" 

OUTPUT 711 3; "ASSSLIMNIMD4SP10AD" ,F1,"SCIAV4NU4AM2" 
J=(V1-A(1))/1I+1 

REM ** DATA ACQUISITION FOR SET DRIVE AMPLITUDE ** 
FOR R=1 TO J 

OUTPUT 711 3;"“RE" 

WAIT 11500 

ENTER 723 3 V 

E(R) +¥V 

REM ** STORE AMPLITUDE OF FUND AND 1ST TWO HARMONICS *+* 
FOR N=1 TO 3 

M= (F (CN) -F1) /1000#250 

OUTPUT 711 3;"MP",M 

OUTPUT 711 3 “LMK" 

ENTER 711 3; B(R,N) ,C(R,N) 

NEXT WN 

D=R+1 

A(D) =A(R) +1 

OUTPUT 7O7 3;"AP",AC(D) , "MV" 

NEXT R 

REM ** OUTPUT DATA TO PRINTER *+* 

PRINTER IS 701,80 

PRINT “FUND FRE@ (HZ) =",F(1) 


PRINT “SOURCE VOLTS LC10" 

PRINT "MILLIVOLTS(P-P) VOLTS(RMS) VF1IC(CRMS) VF2¢(RMS) 
VFS (RMS) “ 

FOR R=1 TO J 


IMAGE 50,12X,D.DDE,4X,D.DDE,4X,D.DDE,4X,D. DDE 
PRINT USING 400 ; A(R) ,E(R),B(R,1),B(R,2) ,B(R,3) 
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470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 


NEXT R 
REM ## RESET INSTRUMENTATION TO PRESET CONDITION *# 
OUTPUT 707 3"PR" 

OUTPUT 711 3"PRS" 

REM ## STORE DATA ON TAPE «* 
ASSIGN# 1 TO "HYD1" 

PRINT# 1,1 3 AC) 

PRINT# 1,2 3 E() 

PRINT# 1,3 3 BC,) 

ASSIGN# 1 TO * 

REM #*# PLOT DATA ON COMPUTER PRINTER ** 
PRINTER IS 2 

FOR K=1 TO 3 

GCLEAR 

SCALE 0,10000,0,.01 

XAXIS 0,1000 

YAXIS 0,.001 

PENUP 

FOR I=1 TO J 

PLOT ACI) ,B(I,K) 

NEXT I 

COPY 

NEXT K 

GCLEAR 

SCALE 0,10000,0,.01 

XAXIS 0,1000 

YAXIS 0,.001 

FOR K=1 TO 3 

PENUP 

FOR I=1 TO J 

PLOT A(I) ,B(I,K) 

NEXT I 

NEXT K 

COPY 

DISP "END" 

END 
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SAMPLE RAW DATA 


PZT #1 

Laser Diode 35 milliamps 

Current 

Fundamental 

frequency (Hz) 72 

Source Fundamental ist Harmonic 2nd Harmonic 
Voltage Voltage Voltage Voltage 
(volts) (millivolts) (millivolts) (millivolts) 
(rms) (rms) (rms) (rms) 
0.035 0.783 0.030 0.0019 
0.033 1.180 0.066 0.0039 
0.071 1.370 0.105 0.0048 
0.071 1.980 0.182 0.0106 
0.089 2.320 0.239 0.0164 
0.106 2.700 0.349 0.0241 
0.124 3.050 0.456 0.0347 
0.142 3.430 0.544 0.0473 
0.1389 3.800 0. 642 0.0647 
0.177 4.070 0.812 0.0830 
0.195 4.380 0.916 0.1150 
0.230 4.820 1.130 0.1460 
0.247 5.010 1.270 0.1720 
0.247 3.340 1.510 0.2060 
0.265 32270 1.610 0.2540 
0.283 5.790 1.750 0.2880 
0.300 6.110 2.010 0.3280 
0.318 6.310 2.150 0.4070 
0.336 6.460 2.370 0.4520 
0.353 6.850 2.440 0.5200 
0.373 6.870 2.460 0.3770 
0.390 72130 2.960 0.6490 
0. 409 7.240 2.970 0.7160 
0.426 7.2300 32 220 0.8160 
0. 444 7.370 3. 420 0. 8880 
0.461 7.2000 3.2630 0.9650 
0.478 7.700 3.830 1.0300 
0.497 7.7950 3.2990 1.1300 
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0.5135 
0.532 
0.549 
0.567 
0.5864 
0.603 
0.620 
0.638 
0.655 
0.672 
0.491 
0.707 
0.726 
0.743 
0.761 
0.778 
0.797 
0.814 
0.832 
0.849 
0.868 
0.885 
0.902 
0.920 
0.937 
0.9356 
0.9735 
0.991 
1.010 
1.030 
1.040 
1.060 
1.080 
1.100 
1.110 
1.1350 
1.1350 
1.170 
1.190 
1.200 
1.220 
1.240 
1.260 
1.270 
1.290 
1.3510 
1.3350 
1.340 
1.3560 
1.380 
1.400 


7.650 
7.750 
7.890 
7.460 
727350 
7.770 
7.450 
72-4350 
7.460 
7.2550 
7.280 
7.070 
6.930 
6.850 
6.590 
6.500 
6.200 
6.050 
43-760 
3.0 
3-2 290 
4-070 
4.770 
4.490 
4.200 
3.870 
3.590 
35-270 
35-000 
2./10 
22560 
2.040 
1.750 
1.480 
1.170 
0.858 
0.617 
0.346 
0.339 
0.611 
0.841 
1.080 
1.390 
1.610 
1.920 
2.170 
2.570 
2.600 
2.8350 
2.9780 
35.260 
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4.070 
4.520 
4.250 
35.990 
4.690 
4.340 
3.030 
4.720 
3-190 
4.910 
~.- 080 
3- 450 
4. 460 
3-650 
v3. 750 
3. 820 
3-970 
6.050 
6.100 
6.520 
6.250 
6.370 
6.500 
6.480 
6.540 
6.540 
&. 480 
6. 480 
6.560 
6.500 
37-910 
6.210 
6.170 
6.120 
6.070 
Je JO 
3-2 6350 
3. 460 
32520 
3-150 
4.840 
4.710 
4.570 
4.280 
4.0350 
3.4350 
3.5960 
3.570 
5.070 
2.840 
2.660 


1.2600 
1.3100 
1.4400 
1.5000 
1.6700 
1.7400 
1.8100 
1.9400 
2.90200 
2.2200 
2.2900 
2.4000 
2.3400 
2.6700 
2.8100 
2.9500 
35.0400 
3.1800 
3.2900 
3.54600 
3.4900 
3.6100 
3.6900 
3.8500 
4.0000 
4.0600 
4.1900 
4.2200 
4.4500 
4.5600 
4.6200 
4.7300 
4.7200 
4.8400 
4.8200 
»- 0100 
3. 1200 
~- 0600 
43. 1600 
4. 2500 
3. 2900 
». 2800 
3.2900 
~- 2800 
3.2900 
3. 5600 
37-4000 
3. 5400 
J. 2900 
3. 2900 
3. 1900 


1.420 
1.430 
1.450 
1.470 
1.490 
1.500 
1.520 
1.540 
1.560 
1.570 
1.590 
1.610 
1.4630 
1.440 
1.460 
1.680 
1.700 
1.720 
1.730 
1.750 
1.770 
1.790 
1.800 
1.820 
1.840 
1.860 
1.870 
1.890 
1.910 
1.930 
1.9350 
1.960 
1.980 
2.000 
2-030 
2.050 
2-970 
2.080 
2.100 
2.120 
2.140 
2.150 
2.170 
2.190 
2.210 
2-210 
222350 
2-240 
2.260 
2.280 
2.500 


3-540 
3.740 
3-420 
3.900 
3 ISO 
4.110 
4.200 
4.230 
4.3640 
4.440 
4.510 
4.560 
4.530 
4.410 
4.560 
4.530 
4.460 
4.3550 
4.250 
4.250 
4.160 
4.010 
3-920 
3.890 
3.780 
3.470 
32950 
3.2210 
3.250 
3.030 
2.820 
2.970 
2.260 
22520 
2-990 
1.940 
1.710 
1.510 
1.270 
1.090 
0.846 
0.4683 
0.474 
0.320 
0.361 
0.524 
0.667 
0.909 
1.080 
1.260 
1.470 
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2.410 
2.210 
1.930 
1.780 
1.500 
1.310 
1.080 
0.863 
0.583 
0.366 
0.293 
0.345 
0.601 
O. 766 
1.020 
1.220 
1.420 
1.580 
1.780 
1.930 
2-090 
22290 
2. 480 
2.720 
2.890 
3.020 
2.880 
323/70 
3.440 
3.630 
3.7350 
3.860 
3.570 
3.990 
3.980 
4.080 
4.010 
4.120 
4.110 
4.050 
4.030 
3.960 
3.830 
35.880 
3.820 
3.760 
3.640 
3.35/70 
3.350 
3-220 
3-090 


~- 1400 
3. O500 
4.8700 
4.9600 
4.8100 
4.7100 
4.5800 
4.4800 
4.3900 
4.2700 
4.1200 
3.9900 
35.8100 
3. 4000 
3.4900 
3.5400 
3.1300 
2.9800 
2.7600 
2.49800 
2.4400 
2.2500 
2.0300 
1.8400 
1.4600 
1.4400 
1.2600 
1.0400 
0.8710 
0.46390 
0.4490 
0.2430 
0.90772 
0.1690 
0.3640 
0.5500 
0. 7450 
0.9230 
1.1200 
1.2800 
1.4200 
1.46300 
1. 7800 
1.9200 
2.9700 
2.1900 
22 3000 
2. 4300 
2.1400 
2.7100 
2.7800 


2.510 
2.3550 
2.2300 
0.370 
2. 580 
2.400 
2.420 
2.440 
2.4640 
2.2470 
2.490 
2.310 
2.2050 
2.240 
2.240 
2.480 
2.600 
2.410 
2.630 
0.4650 
2.4/0 
2.680 
2.790 
2.720 
2.740 
2.740 
2.7709 
2.790 
2.810 
2.830 
2.840 
2.860 
2.880 
2.900 
2.920 
2.9730 
2.9750 
2.970 
2.2990 
3.900 
3.920 
3.040 
3.960 
3.970 
3.990 
3.2110 
3.2130 
3.2140 
3.140 
5.180 
3.200 


1.440 
1.830 
1.990 
1.970 
2.190 
2.430 
2.970 
2.790 
2.790 
2.900 
0.950 
3.130 
3.180 
3.180 
322/09 
3.2430 
3.580 
3-430 
3.430 
3.450 
3.430 
3.23590 
3.490 
3.470 
32590 
3.320 
3.2310 
3.190 
3.1350 
3.050 
2.970 
2.870 
2.740 
2.420 
2.419 
2.2590 
2.2/0 
2.1350 
1.990 
1.790 
1.490 
1.520 
1.360 
1.200 
0.997 
0.855 
0. 686 
0.558 
0. 404 
0.303 
0. 454 
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2.9740 
2.810 
2.490 
2.920 
2.25/70 
2.160 
2.010 
1.870 
1.740 
1.540 
1.390 
1.210 
1.030 
0.746 
0.678 
0.452 
0.320 
0.259 
0.253 
0.292 
0. 438 
0.586 
0. 800 
0.939 
1.140 
1.240 
1.510 
1.420 
1.820 
1.970 
2.130 
2.260 
2.3500 
2.740 
2.400 
2.720 
2.780 
2.9710 
3.000 
3.020 
3.100 
3.970 
3.070 
3.250 
3.250 
3.220 
3.220 
3.2130 
3.100 
3.970 
3.020 


2.7000 
3.0400 
3.0900 
3.20700 
3.1900 
3.2100 
3.3000 
3.5200 
3.3600 
3.353400 
3.2900 
3.3700 
3.35400 
3.3100 
3.23100 
3.3000 
3.1600 
3.1000 
3.0500 
3.90500 
22-9700 
2.8700 
2.77900 
2. 4800 
2.3200 
2.4300 
2.5200 
2.2100 
2.0500 
1.9100 
1. 7800 
1.4500 
1. 4800 
1.35400 
1.1900 
1.0500 
0.8880 
0.7160 
0.5690 
0.3920 
0.2440 
0.0733 
0.0791 
0.2470 
0. 3600 
0.5340 
0.6920 
0.8300 
0.9820 
1.1300 
1.2800 


3.210 
3-230 
3. 250 
3-270 
3.290 
3.3500 
3-520 
3-340 
3-2 350 
3.2370 
32 390 
3.410 
3.430 
3.440 
3-460 
3.480 
3.300 
3.020 


0.585 
0.748 
0.904 
1.030 
1.210 
1.330 
1.480 
1.550 
1.730 
1.870 
1.990 
2.130 
2.270 
2.350 
2.490 
2-470 
2.260 
2.650 
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3.000 
2.720 
2.860 
2.640 
2.680 
225920 
2-420 
223350 
2.150 
2.970 
1.940 
1.850 
1.710 
1.560 
1.410 
1.280 
1.130 
0.984 


1.4000 
1.5200 
1.6400 
1.7800 
1.8700 
1.9500 
2.0400 
2.1500 
2.1800 
2.5100 
2.5500 
2.4500 
2.9500 
2.4400 
2.6400 
2.46000 
2.6500 
2.6300 
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